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ABSTRACT 

Point  mass  calculations  are  used  to  model  switched  implosions  on  several  pulsed  power  machines.  The  model 
includes  a  lumped  circuit  representation  of  the  pulsed  power  source.  A  simple  switching  model  is  used  to  describe  a 
standard  plasma  flow  switch.  Implosion  kinetic  energies  are  obtained  at  a  convergence  ratio  of  20  to  1.  Heuristic 
arguments  are  used  to  estimate  the  plasma  temperature  at  pinch,  the  total  x-ray  output  and  the  radiation  pulse  width. 
Switched  models  are  presented  for  Pegasus  II,  Shiva  Star,  Procyon  and  Atlas. 

I.  INTRODUCTION 

A  simple  point  mass  model  is  used  to  calculate  switched  implosions  on  a  variety  of  pulsed  power  machines.  The 
kinetic  energy  and  velocity  of  the  load  at  a  convergence  of  20: 1  are  used  to  estimate  die  plasma  temperature  of  a 
pinch,  die  total  x-ray  output  and  the  pulse  width.  The  discussion  includes  proposed  Atlas  parameters,  but  results  for 
Pegasus  II  and  Shiva  Star  switched  implosions,  as  well  as  for  Procyon  high  explosive  (HE)  driven  generator  driven 
implosions  are  included.  For  some  of  these  experimental  data  and  extensive  2-D  models  exist  and  may  be  used  to 
put  die  Atlas  results  in  perspective.  The  point  mass  model  couples  a  single  loop,  lumped  external  circuit  equation 
for  a  capacitor  bank  or  HE  generator  to  dynamic  equations  describing  the  load  implosion.  A  dynamic  model  is 
included  for  die  motion  of  die  switch  down  die  coax  and  for  the  transfer  of  current  to  the  load  (switching).  The 
equations  are  expressed  in  units  of  g,  cm  and  ps.  In  diis  system  die  natural  units  include:  current  in  10^  amps 
(DMA);  voltage  in  10^  volts  (DkV);  inductance  in  nH;  capacitance  in  mF;  and  resistance  in  m£2.  The  unit  of 
energy  is  die  energy  unit,  1  eu  =  10^  ergs  =  0.1  M.L  Computational  results  will  be  expressed  in  conventional 
units.  The  pressure  gradient  in  the  plasma  which  produces  bounce  is  not  included  in  die  point  mass  model.  * 
Therefore,  die  calculations  are  followed  only  to  a  radius  of  about  0.25  cm,  corresponding  to  a  convergence  of  20:1. 
This  gives  a  pinch  radius  typical  of  some  two-dimensional  perturbed  implosion  models. 
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Figure  1.  Single  loop  equivalent  circuit  model 
for  capacidvely  driven  pulsed  power  system. 
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Figure  2.  Single  loop  equivalent  circuit  model 
HE  driven  pulsed  power  system  (Procyon). 


Specific  inducdve  store  systems  are  represented  by  die  equivalent  circuits  shown  in  Figure  1  (Atlas  parameters 
for  die  capacitively  driven  system)  and  Figure  2  (Procyon  parameters  for  the  HE  system).  System  specific  values 
for  die  circuit  parameters  will  be  defined  below.  The  bank  or  HE  system  is  modeled  by  a  single  loop  potential  drop 
equadon  which  relates  die  bank  parameters  or  storage  inductor  parameters  to  die  current  and  the  change  in  inductance 
of  die  load.  The  inductance  change  of  the  switch  is  also  included.  The  damping  resistor  has  not  been  included, 


t  Permanent  address:  Department  of  Chemical  and  Nuclear  Engineering,  University  of  New  Mexico,  Albuquerque, 
New  Mexico. 

^Simple  terms  can  be  included  in  die  model  that  would  produce  a  bounce  and  pinch.  The  most  physical  is  associated 
widi  a  axial  magnetic  field  seeded  inside  the  load  foil  at  die  beginning  of  die  implosion.  Odier,  lieurisdc  terms  can 
also  be  constructed  to  represent  die  plasma  pressure.  These  do  not  include  die  effects  of  shocks  formed  when  die 
plasma  stagnates  on  axis,  and  thus  should  not  be  considered  reliable. 
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although  the  equations  below  could  be  easily  modified  to  include  it.  The  current  predicted  without  the  damping 
resistor  will  therefore  be  an  overestimate. 

The  point  mass  model  discussed  below  is  useful  in  understanding  the  coupling  between  prime  power  source, 
switch  and  load  and  the  overall  performance  characteristics  of  the  implosions.  Extensive  numerical  models  of  these 
phenomena  using  1-D  and  2-D  codes,  as  well  as  data  from  different  pulsed  power  regimes,  has  shown  that  the 
behavior  of  the  implosion  is  strongly  dominated  by  instability  effects.  We  have  used  the  results  of  these  studies  as  a 
guide  in  picking,  for  example,  the  convergence  ratio  at  which  to  evaluate  the  radiation  output  from  the  model.  In 
addition  to  their  heuristic  value,  a  comparison  of  the  models  below  may  be  useful  in  studying  trends. 

Section  II  describes  the  point  mass  model  for  switched  implosions  in  general.  Section  III  discusses  the  results 
from  switched  implosions  on  Pegasus  II,  Shiva  Star,  Procyon  and  Atlas.  Section  IV  summarizes  the  results. 


II.  POINT  MASS  MODEL 

The  point  mass  ntodel^  consists  of  a  lumped  circuit  equation  for  the  pulsed  power  system  and  the  dynamic  equations 
for  the  switch  and  the  load.  We  consider  a  capacitively  driven  system  first.  The  circuit  equation  is 

Q  ^  d(LI) 

C  dt 


+  IR  =  0 


(l) 


where  Q  is  the  instantaneous  charge  on  the  bank,  C  is  its  capacitance,  I  is  the  current,  R  is  the  bank  resistance  and  L 
is  the  total  inductance  in  the  system.  The  latter  is  given  by 


L(t)  -  L0  +  LPFS(t)  +  f(z)  LlvaJt) 


(2) 
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Here  L0  is  the  inductance  of 
the  bank  and  of  the  power 
flow  channel,  and  LpFS  is 
the  inductance  associated 
with  the  instantaneous 
position  of  the  switch 
plasma.  The  last  term  on 
the  right  side  of  (2)  contains 
two  quantities:  the  total 
inductance  of  the  load 
region,  Lioad,  which  is 
defined  by  the  instantaneous 
position  of  the  load;  and  a 
switching  function,  f(z), 
which  describes  the 

„.  oo.  .  ,  .  ,  „  .  ,  ,  ...  nr,c  instantaneous  fraction  of  tlie 

Figure  3.  Schematic  ol  a  typical  power  flow  channel  and  coax  containing  a  PFS  bank  currem  that  has  been 

and  a  load  foil.  The  PFS  is  located  initially  at  z0  and  moves  to  the  left.  transferred  to  the  load  The 

product  describes  the 

inductive  change  associated  with  switching  into  the  load  region.  The  inductance  Lpps  depends  on  the  switch's 
motion,  and  is  given  by  (see  Figure  2) 


-PFS 


=  2  (z0-z(t))  ln(~\ 

\  / 


(3) 


were  zq  is  the  initial  position  of  the  PFS,  z(t)  is  its  instantaneous  position  in  the  coax,  and  Rj  and  Rq  are  the  inner 
and  the  outer  radii  of  the  coax,  respectively.  In  defining  (3),  the  motion  of  the  coax  is  understood  to  be  from  right  to 
left  with  respect  to  its  original  position.  For  all  of  the  examples  below  Rj  =  7.62  cm  and  Ro  =  10.16  cm. 


II.A.  Switching 

Two-dimensional  models  of  a  standard  PFS  indicates  approximately  planar  motion  of  the  switch  plasma  down  the 
coax.  When  the  switch  moves  off  the  up-stream  edge  of  the  load  region,  magnetic  flux  is  transferred  to  the  load  on 


^We  assume  a  "standard"  PFS,  as  described,  for  example,  in  P.J.  Turchi,  M.L.  Alme,  F.  Bird,  C.N.  Boyer,  S.K. 
Coffey,  D.  Conte,  J.F.  Davis  III,  and  S.W.  Seiler,  IEEE  Transactions  on  Plasma  Science,  Vol.  PS-15,  No.  6,  1987, 
pg.  747. 
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essentially  the  Alfven  time  scaled  Assuming  that  the  switch  remains  open  during  the  implosion,  the  only  effect  of 
the  PFS  after  switching  is  to  produce  a  slight  additional  inductance  as  it  moves  across  the  load  slot  and  the  field  is 
stripped  off  by  the  down-stream  vanes  (see  Figure  3).  Thus  switching  may  be  modeled  in  die  following  way.  First, 
when  die  PFS  open  into  the  load  region  diere  is  an  increase  in  inductance  which  is  associated  with  die  volume 
between  die  PFS  and  die  initial  position  of  die  load.  This  will  occur  over  some  time  scale,  which  is  input  for  diis 
model.  We  assume  diat  die  inductance  change  is  given  by  the  volume  bounded  by  die  inner  radius  of  die  electrode 
and  die  instantaneous  position  of  the  load: 

Lbad  =  2  Az  (4) 

Here  Az  is  die  widdi  of  the  load  slot  (Figure  2),  and  r(0)  is  die  initial  position  of  die  load  foil  To  model  switching 
we  multiply  diis  inductance  by  a  form  factor,  f(z(t)),  which  depends  on  die  instantaneous  location  of  die  PFS,  z(t), 
and  on  a  lengdi  scaled  Azs  over  which  complete  transfer  of  current  to  die  load  region  occurs: 

f(z(f))  =  %  Azj  ~  (5) 

Here  %  is  an  efficiency  factor”  (unity  for  complete  switching).  The  parameter  z\  represents  die  axial  point  at  which 
die  switch  begins  to  open.  Typically  this  will  be  taken  as  the  up-stream  edge  of  die  load  region.  Thus  f(z)  is 
essentially  zero  for  z  far  from  z\,  and  dien  transitions  to  unity  as  die  PFS  moves  a  characteristic  distance  Azs.  The 
form  of  die  switching  function  (5)  has  been  chosen  for  convenience.  Odier  functions  which  vary  smoodily  from 
zero  to  one  could  also  be  used.7  The  switching  time  is  dien  determined  by  diis  distance  and  die  velocity  of  die  PFS 
as  it  crosses  die  load  slot.  Since  the  load  does  not  move  significantly  during  diis  time,^  die  radius  r(t)  in  (4)  is 
essentially  constant  until  switching  is  nearly  complete^.  Note  diat  die  time  derivative  of  f(z)  is  non-zero  during 
switching,  so  diat  dLioa(i  /  dt  is  non-zero  even  if  the  load  radius  r  is  a  constant. 

As  die  PFS  moves  down  die  coax  die  inductance  behind  it  increases,  and  diis  reduces  die  current  diat  can  be 
supplied  by  die  bank.  We  assume  diat  the  bank  inductance  corresponds  to  all  of  die  volume  up  to  die  back  of  die 
PFS  .  When  switching  begins  diere  will  be  a  furdier  decrease  in  current.  Finally,  as  die  implosion  proceeds,  die 
load  inductance  will  increase  as  die  radius  of  die  load  decreases,  and  a  furdier  decrease  in  current  will  occur. 

Once  die  bank  current  has  been  obtained  die  load  current  is  obtained  as  die  product  of  die  switching  function  and 
die  bank  current: 

IL(t)  =  f(z(t))  I(t)  (6) 

Thus  die  load  current  follows  from  a  consistent  solution  of  die  circuit  equation  and  die  dynamic  equation  for  die 
switch  mass.  Note  that  die  latter  also  depends  on  die  bank  current.  As  the  bank  current  decreases,  so  will  die  load 
current. 


II.  11.  Circuit  Equation 

Substituting  die  results  above  into  the  circuit  equation  (1),  and  expressing  die  current  in  terms  of  die  charge  on  die 
bank,  gives  die  equation  for  die  circuit  coupled  to  die  inductive  changes  associated  with  die  PFS  and  die  load: 

2 

0  =  ®  +  {Lc  +  LPFS(t)  +  f(z)  Lloa^—Q-  + 


^ This  aspect  of  die  point  mass  switching  model  is  contained  in  die  parameters  appearing  in  f(z)  defined  in  (5)  below. 

^In  die  model  r(0)  =  r(ts),  where  ts  is  die  time  at  which  switching  begins. 

^This  should  not  be  confused  widi  die  height  of  die  load  slot,  which  is  denoted  by  Az 

^Typically  we  expect  0.9  <  q  <  1.0.  We  also  assume  diat  once  die  switch  has  opened  it  remains  open.  This  aspect 
of  die  model  could  be  relaxed. 

7A  physically  based  model  for  the  switch  function,  based  for  example  on  a  time  varying  resistivity,  could  be  developed  for  use 
here. 

^Vaporization  of  die  load  foil  occurs  during  diis  switching  period,  and  actual  inward  radial  motion  tends  to  occur 
only  after  most  of  die  current  has  been  transferred  to  die  load. 

^The  model  actually  does  allow  die  load  to  move  during  switching.  Typical  calculations  show  that  for  load  masses 
of  interest  very  little  motion  begins  before  switching  is  complete. 


1150 


This  represents  the  complete  circuit  equation.  It  depends  on  the  dependent  variables  z(t),  describing  the  motion  of  the 
PFS  down  die  coax,  the  switch  function  f(z),  and  r(t)  describing  die  motion  of  die  load  during  die  implosion.  The 
resistance  of  die  system  could  be  extended  to  include  time  dependent  effects,  such  as  a  fuse  or  die  variable  resistance 
of  die  load  during  die  implosion.  For  die  systems  considered  below  die  load  resistance  is  approximately  R  =  t|l  /  A, 
where  1  =  2  cm  is  die  load  height,  and  A  =  2?trArs  is  die  cross  secdon  area  of  die  current  carrying  region.  When  die 
resishvity  t|  is  largest,  A  ~  5  cm^  and  R  is  less  dian  about  0.05  mfl.  The  actual  value  varies  during  die  implosion, 
but  is  generally  small  compared  widi  die  bank  resistance.  Dynamic  equations  for  die  variables  appearing  in  (7)  are 
described  next. 


II. C.  Equations  of  Motion 

For  planar  modon  of  an  ideal  PFS  down  the  coax,  die  z  component  of  die  acceleradon  will  be  independent  of  radius 
and  die  radial  acceleration  will  vanish.  Using  the  relation  for  a  toroidal  magnetic  field  2  I  =  r  Bq,  die  point  mass 
axial  equation  of  motion  is 


mPFS 


d2z 

dt2 


(8) 


where  nipps  is  die  total  mass  of  the  PFS.  This  assumes  that  die  mass  distribution  is  close  to  Mr2-,  so  diat  die 
axial  motion  is  essentially  planar.  The  dependent  variables  are  die  current,  I(t),  and  die  axial  position,  z(t).  No 
provision  has  been  made  to  model  flux  stripping  as  die  plasma  moves  past  die  vanes  and  out  of  die  coax.  10  in  tJic 
models  below  die  PFS  usually  has  not  moved  far  enough  during  die  implosion  of  die  load  for  tliis  to  be  necessary. 

The  radial  equation  of  motion  for  the  load,  whose  mass  is  denoted  by  nip,  is 


mL 


£r 

dt2 


,2  Az 

lL~ 


(9) 


where  Az  is  die  height  of  die  load  foil,  and  Il(0  is  the  current  in  die  load.  The  dependent  variables  are  die  radius  r(t), 
and  die  current  Il(0.  As  noted  previously,  the  model  does  not  include  a  plasma  pressure  gradient  which  would  be 
needed  to  simulate  the  bounce  at  pinch. 


II. D.  High  Explosive  Generator  Driven  Systems 
The  Procyon  system,  which  has  been  described  elsewhere  1 1 ,  is  modeled  by  a  single  loop  circuit  containing  a  storage 
inductance,  Lst,  a  time  dependent  closing  switch  S(t),  and  a  power  flow  inductance  L0  (Figure  1).  The  models  above 
are  used  for  die  PFS  and  the  load.  An  initial  current,  I0,  is  set  in  the  storage  inductor  and  die  switch  S  is  open.  The 
problem  starts  by  ramping  up  a  closing  switch  function  S(t),  which  delivers  current  to  die  PFS  widi  a  10%  to  90% 
rise-time  of  about  1.6  ps  as  observed  in  the  PC3  experiment.  The  closing  switch  function,  fs(t)  is  chosen  to  be  of 
die  form 


fs(f)  -  e[(t-tJ/ A'l  +  j) 


(10) 


The  switch  begins  to  close  at  to,  which  we  set  to  zero,  and  At  represents  the  time  scale  for  closing.  The  factor 
represents  efficiency  of  closing.  For  simplicity  we  set  the  resistance  of  die  system,  R  =  0.  The  equivalent  circuit 
equation  is  dien 


L  dl 

s,dt 


+ 


d(Ut)l) 

dt 


0 


(11) 


For  die  HE  system,  L(t)  includes  the  inductance  of  the  PFS,  die  power  flow  channel,  die  load,  and  die  effects  of  die 
closing  switch  S(t).  The  current  in  die  PFS  is  given  by  fs(t)  I,  and  die  current  in  die  load  by  fs(t)  f(z)  I.  In 
evaluating  (1 1)  we  include  die  rate  of  change  of  bodi  fs(t)  and  f(z). 

Finally,  we  note  diat  the  models  above  can  be  easily  modified  to  describe  capacitively  driven  or  HE  driven 


l^A  simple  stripping  model  could  be  included  using  a  form  factor,  fs(z),  which  would  supply  a  right-ward  acting 
force  on  die  PFS  starting  at  a  specified  point  in  the  coax,  and  bringing  the  plasma  to  a  halt  over  some  lengdi  scale. 
Such  a  term  would  be  included  on  die  right  side  of  (9)  to  reduce  the  axial  acceleration. 

^J.  Gofordi,  et.  al„  Megagauss  Magnetic  Field.  Generation  and  Pulsed  Power  Applications,  Ed.  by  M.  Cohen  and 
R.B.  Spielntan,  (Nova  Science  Publishers:  New  York;  1994),  pg.  841;  J.  Gofordi  et.  al.  Proceedings  of  the  Ninth 
IEEE  International  Pulsed  Power  Conference,  Albuquerque,  June  1993,  pg.  36;  R.L.  Bowers,  A.E.  Greene,  D.L 
Peterson  and  N.F.  Roderick,  Proceedings  of  the  Ninth  IEEE  International  Pulsed  Power  Conference,  Albuquerque, 
June  1993,  pg.  538. 
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include:  interaction  of  the  switch  plasma  and  the  electrodes  as  it  moves  down  the  coaxial  barrel;  the  effects  of 
distributed,  background  plasma  density  on  die  switch's  efficiency  and  its  ability  to  remain  open  during  the  implosion 
and  pinch;  effects  of  magnetically  driven  Rayleigli-Taylor  and  wall  induced  instabilities  on  die  load  plasma; 
instabilities  induced  by  die  switch  plasma  during  switching;  and  absorption  of  pinch  radiauon  by  switch  plasma. 

Table  1.  Input  Parameters  for  Switched  Implosion  Models 


Regime 

Voa 

(kV) 

Cb 

(mF) 

zi 

(cm) 

Az 

(cm) 

rnpps 

(mg) 

mL 

(mg) 

R 

(mn) 

LPFC 

(nH) 

Azs 

(cm) 

Pegasus  II 

84 

0.85 

4.45 

4.45 

60 

12 

0.5 

23 

0.95 

Shiva  Star 

90c 

1.313 

6.3 

6.3 

120 

20 

1.0 

19 

0.10 

0.95 

90d 

1.313 

6.3 

6.3 

120 

12.7 

1.0 

19 

0.28 

0.95 

Procyon 

21.7 

79 

3.5 

6.3 

200 

30 

0. 

28 

0.30 

0.95 

Adas 

540 

0.2 

6.3 

6.3 

200 

14 

8.0 

23 

0.05 

0.95 

For  Procyon  a)  is  initial  current  (MA)  in  storage  inductor;  b)  storage  inductance  (nH);  c)  20  mg 
load;  d)  13  mg  load. 

Table  2.  Comparison  of  Switched  Implosion  results. 


Regime 

Erada 

(MJ) 

dE/dt a 
(MJ/ps) 

Tcb 

(eV) 

Imax 

(MA) 

x  c 
(ns) 

*s 

(ns) 

Teff 

(eV) 

Pegasus  II  (12  mg) 

0.45 

22 

190  -  320 

9.2 

20 

50 

116 

0.35 

18 

145  -  250 

9.2 

20 

200 

Shiva  Star  (20  mg) 

0.63 

16 

170  -  280 

12.2 

39 

70 

(13  mg) 

0.67 

22 

260  -  440 

12.2 

31 

70 

116 

Procyon  (30  mg) 

1.2 

67 

200  -  330 

13.8 

18 

350 

152 

Atlas  (14  mg) 

4.5 

750 

1600  -  2700 

32.6 

6 

20 

280 

a)  Assumes  50  %  of  implosion  KE  goes  into  radiation;  b)  Uses  (12),  with  Cy  =  1200  -  2000  eu  g'^keV’^; 


c)  Calculated  as  twice  0.25  cm  /  vr,  where  vr  is  die  load  velocity  at  0.25  cm. 


Figure  4e.  Current  in  die  PFS  and  in  die  load  versus 
dme.  The  onset  of  switching  is  shown  is  shown  by 
die  switch  function  f(t)  in  Figure  4d  above. 


Figure  4f.  Implosion  kinetic  energy  in  MJ,  versus 
time.  The  kinetic  energy  at  20: 1  convergence  is 
about  4.5  MJ. 


Applications,  Ed.  by  M.  Cohen  and  R.B.  Spiehnan,  (Nova  Science  Publishers:  New  York;  1994),  pg.  713. 
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direct  drive  experiments. 


III.  OUTPUT  FROM  SWITCHED  IMPLOSIONS 

Although  the  point  mass  model  does  not  include  radiative  losses,  we  can  make  order  of  magnitude  estimates 
for  quantities  of  interest  from  the  results  above  and  several  simple  assumptions^.  Starting  with  die  relationship 
between  die  specific  kineuc  energy  of  die  imploding  load  and  its  diennal  energy  per  gram 

CVT  =  (1/2)  v2  (12) 

we  can  estimate  die  central  temperature  of  die  pinch.  Here  Cy  is  a  pseudo-heat  capacity,  whose  value  is  known 
from  detailed  calculations.  We  use  values  for  die  pseudo-heat  capacity  in  die  range  1200  -  2000  eu  g'l  keV'l  (10^  - 
1.7  x  10s  erg  g'^  K'  1 ),  and  die  load  velocity  at  20:1  convergence.  Results  of  two-dimensional  calculations  indicate 
that  die  amount  of  energy  radiated  often  equals  or  exceeds  die  kinetic  energy  near  pinch.  We  shall  assume  that  all  of 
die  kinetic  energy  at  20: 1  convergence  goes  into  radiation,  and  set  q  =  1. 

As  a  specific  example,  consider  an  Atlas  switched  implosion  of  a  2  cm  high,  5  cm  radius  load  widi  a  mass  of 
14  mg.  We  assume  a  200  mg  switch  plasma,  initially  located  6.3  cm  to  the  right  of  die  upstream  edge  of  die  load 
slot  (z0  -  z\  in  Figure  2).  The  switch  mass  was  arrived  at  by  requiring  diat  die  PFS  reach  die  load  region  witii  an 
axial  velocity  of  8  -  10  cm  /  ps,  but  before  peak  current  delivery  into  a  dead  short.  The  initial  bank  voltage  was  V0 
=  540  kV  and  C  =  0.2  mF,  corresponding  to  a  bank  energy  of  29.2  MJ.  The  inductance  of  die  bank  and  die  power 
flow  channel  up  to  die  back  of  the  initial  position  of  die  PFS  is  assumed  to  be  23  tiH.  The  switch  moves  to  die  left 
toward  die  load  region  as  in  Figure  2.  The  onset  of  die  pinch  is  taken  to  be  a  load  radius  of  0.25  cm.  The  switch 
function  (5)  uses  Azs  =  0.05  cm,  £  =  0.95  for  95%  switching  efficiency,  zq  =  0.0  cm,  and  zj  =  -6.3  cm.  The  point 
mass  model  predicts  an  implosion  time  to  a  20:1  convergence  of  2.017  ps.  The  current  flowing  in  die  PFS  reaches 
a  peak  value  of  32.3  MA  prior  to  switching  into  die  load  region.  Switching  starts  at  ts  =  1.82  ps,  and  is  complete  in 
about  40  ns.  The  10%-90%  current  rise  time  in  die  load  is  about  20  ns. 

The  actual  implosion  time,  as  measured  to  a  radius  of  0.25  cm,  is  about  160  ns.  The  load  velocity  reaches 
about  80  cm/ps  at  0.25  cm,  which  corresponds  to  a  kinetic  energy  of  4.5  MJ.  A  stagnation  time  for  die  pinch  may 
be  defined  as  2x,  where  x  is  0.25  cm  divided  by  die  load  velocity  at  that  point.  For  Atlas  x  =  2  x  0.25  cm  /  80 
cm/ps  =  .006  ps.  The  total  power  associated  with  die  implosion  is  q  dE/dt  =  4.5  q  MJ  /  0.006  ps  =  750  q  MJ/ps. 
Assuming  a  conversion  efficiency  q  =  1  we  have  Era(j  =  4.5  MJ  and  dE/dt  =  750  MJ/ps.  The  x-ray  output  from  die 
switched  implosion  represents  about  15  %  of  the  initial  bank  energy.  Using  (12)  we  can  estimate  die  temperature  of 
die  pinch  to  be  in  die  range  1.6  -  2.7  keV.  In  this  model  the  peak  current  in  die  load  is  30.6  MA.  A  summary  of 
results  from  die  point  mass  model  is  given  in  Figure  4. 

The  model  has  been  used  to  evaluate  the  performance  of  odier  pulsed  power  systems.  The  parameters  for  each 
are  listed  in  Table  1 .  Table  2  compares  output  and  power  for  Atlas  switched  implosions  widi  results  for  odier  pulsed 
power  systems  in  Table  1 .  It  also  lists  estimates  of  die  pinch  temperature,  Tc,  peak  current  in  die  load,  Imax» 
stagnation  time,  Xc,  switching  time,  xs,  and  die  surface  temperature  of  die  pinch.  The  latter  is  defined  by 

T'S  ■  30  *v 

where  die  surface  radius  reff  =  0.25  cm,  a  is  the  Stefan-Boltzmann  constant  and  E  is  in  eu. 

IV.  CONCLUSIONS 

The  point  mass  model  for  switched  implosions  has  been  applied  to  pulsed  power  systems  whose  stored  energy 
ranges  from  3.0  MJ  to  29.2  MJ.  For  most  of  these  systems  2-D  calculations  are  available,  and  have  been  used  to 
motivate  die  choice  of  parameters  for  the  point  mass  models.  The  total  radiation  output  predicted  by  die  model  is  in 
reasonable  agreement  with  experiments  or  calculations.  The  diermalization  times  are  about  10  times  too  short, 
which  results  in  estimated  powers  that  are  large  by  about  die  same  factor.  The  absolute  value  of  die  plasma 
temperatures  are  also  expected  to  be  too  high. 

Two-dimensional  radiation  MHD  models  of  switched  implosions  identify  several  features  which  can  have  a 
dramatic  effect  on  die  quality  of  a  switched  implosion  tiiat  are  not  included  in  die  point  mass  model.  ^  These 


l^These  are  based  in  part  on  experience  with  2-D  perturbed  models  of  implosions.  We  cannot  place  too  much 
credence  in  die  quantitative  nature  of  diese  estimates,  but  die  trends  exhibited  between  die  systems  modeled  below 
may  be  useful. 

^F.J.  Wysocki,  et.  al.,  Megagauss  Magnetic  Field  Generation  and  Pulsed  Power  Applications,  Ed.  by  M.  Cohen 
and  R.B.  Spielman,  (Nova  Science  Publishers:  New  York;  1994),  pg.  767;  R.L.  Bowers,  J.H.  Brownell,  A.E. 
Greene,  D.L.  Peterson,  N.F.  Roderick  and  P.  Turclii,  Megagauss  Magnetic  Field  Generation  and  Pulsed  Power 
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